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ABSTRACT 

Understanding the nature of distant Lya nebulae, oka "blobs," and connecting them to their present- 
day descendants requires constraining their number density, clustering, and large-scale environment. 
To measure these basic quantities, we conduct a deep narrowband imaging survey in four different 
fields, Chandra Deep Field South (CDFS), Chandra Deep Field North (CDFN), and two COSMOS 
subfields, for a total survey area of 1.2 deg 2 . We discover 25 blobs at z = 2.3 with Lya luminosities of 
Lhya= 0.7-8 x 10 43 ergss -1 and isophotal areas of A iso = 10 - 60 U" . The transition from compact 
Lya emitters (A- lso ~a few U") to extended Lya blobs (A- lso > 10 U") is continuous, suggesting a 
single family perhaps governed by similar emission mechanisms. Surprisingly, most blobs (16/25) are 
in one survey field, the CDFS. The six brightest, largest blobs with L\, y(x > 1.5xl0 43 ergss -1 and 
Ai SO > 16 lie only in the CDFS. These large, bright blobs have a field-to-field variance of a v > 1.5 
(150%) about their number density n ~ I.O^qJx 10 -5 Mpc -3 . This variance is large, significantly 
higher than that of unresolved Lya emitters (<j v ~ 0.3 or 30%), and can adversely affect comparisons 
of blob number densities and luminosity functions among different surveys. Our deep, blind survey 
allows us to construct a reliable blob luminosity function. We compare the statistics of our blobs with 
dark matter halos in a 1 h~ x Gpc cosmological N-body simulation. At z = 2.3, n implies that each 
bright, large blob could occupy a halo of M na io ^ 10 13 M if most halos have detectable blobs. The 
predicted variance in n is consistent with that observed and corresponds to a bias of ~7. Blob halos 
lie at the high end of the halo mass distribution at z = 2.3 and are likely to evolve into the ~10 14 M 
halos typical of galaxy clusters today. On larger scales of ~ 10 co-moving Mpc, blobs cluster where 
compact Lya emitters do, indicating that blobs lie in coherent, highly over dense structures. 
Subject headings: galaxies: formation — galaxies: high-redshift — intergalactic medium 
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1. INTRODUCTION 

Lya nebulae, or "blobs," are extended 
z ~ 2-6 with typical Lya sizes of >5 A 
and line lumino sities of Z/Lvr* ^ 10 43 ergss -1 (e.g., 
Keel et al.lll999[ ISteidel et all 12000: iFrancis et al.l [2 001: 
Matsuda et al.l 120041: iDev et al.l 120051: ISmith fc Jarvisl 
20071: iH ennawi et al.l I2009E iPrescott et all 120091 : 
Yang et all [2009). Because the large spatial extent 



of their Lya-emitting gas implies an interaction between 
the surrounding intergalactic medium and any em- 
bedded galaxies, blobs may signal an important phase 
of galaxy formatio n in the early universe, inc l uding 
cold gas acc retion ( Haiman. Spaans. fc Quataert 2000; 
Fard al et al.l 120011 : lYang et al.l "E)Q6: Diikstr a fc Loebl 
2009 ), galactic-scale feed back due to stellar winds 
(jTaniguchi fc Shiova[j2000h . or intense radiative feed- 
back from AGN ([Geach et al.l l2QQ9h . Despite the 
importance of blobs and the controversy regarding their 
origins, even basic properties such as their number 
density, clustering, and large-scale environment are 
poorly constrained. 

To understand into what these mysterious objects will 
evolve in the present day universe, measuring their statis- 
tics is critical due to the direct connection of number 
density and field-to-field variance to halo mass in ACDM 
cosmology. Currently, the halo mass of blobs is unknown. 
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Using the spat i al ext ent and line-width of the Lya line, 
IMatsuda et a l. (2006) estimate dynamical masses of 0.5 - 
20 x 10 12 Mq if the extended Lya emission is from grav- 
itationally bound gas clouds and the resonant scattering 
of Lya can be ignored. 

The clustering of blobs has not been measured directly 
in past work, yet there are hints that it is strong. Af- 
ter surveying over ^4.8 deg 2 in th e NOAO Dee p - Wide 
Bootes field (jJannuzi fc Devlll999D . lYang et al.l (|2009f ) 
discover just four bright blobs, yet two of them lie within 
only ~70" of each other (see also the discovery of a 
Lya blob near a radio-loud Lya h alo by [M atsuda et al.l 
2009) . Some, bu t not all (e.g., IGronwall etHI 120071 : 



Nilsson et al.l 120091 ), narrowband surveys targeting com- 
pact Lya emitters also detect blobs at similar red- 
shifts. For exa mple, while followin g up the two bright 
blobs found by ISteidel et all (|2000h in the SSA22 field, 
IMatsuda et al.l ( 2004f ) discover a spectacular clustering 
of 35 blobs, in which the two brightest lie at the in- 
terse ction of filaments tra ce d by compact Lya emit- 
ters (IMatsuda et al.l I2005D . iPalunas et all (|2004fl dis- 
cover twcQ additional Lya blobs in the J2 143-4423 re- 
gion defined by an o verdensity of compact Lya emitters 
(Fr ancis eFaIll200TI ). To measure directly the number 
density and clustering of Lya blobs, and thus to con- 
strain their halo mass, requires a large volume survey, 
particularly over different sight lines, to account for any 
field-to-field variations. Furthermore, one must apply 

4 One of these blob candidates has now been identified as a low-z 
interloper (J. Colbert priv. comm.) 
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uniform selection criteria for identifying blobs over the 
entire survey volume. 

To acquire a large, unbiased sample of Lya blobs at 
z = 2.3, we have pursued two complimentary narrow- 
band imaging surveys. The shallow, but wide sky 
coverage, survey using the Steward Observatory Bok 
2.3m + 90Prime imager targets r are, luminous Lya 
blobs (L Lya > 2 x 10 43 ergss- 1 ; lYang et al.l I2QQ9I ). 
In this paper, we report the first results from our 
deeper, but smaller sky coverage, survey with the NOAO 
4m telescopes and MOSAIC imagers that targets pre- 
sumably more common intermediate siz e and luminos- 
ity L ya blobs like those discovered by Mats uda et al.l 
(|2004f ). Using blob statistics from the four differ- 
ent 30' x 30' s urvey fields, in the Chandra Deep Field- 
South (CDFS; Brandt et al. 2001), Chandra Deep Field- 
North (CDFN; iGiacconi et al.ll2002[ h and two regions of 
the Cosmic Evolution Survey (COSMOS; iScoville et al.l 
120071 : iKoekemoer et al-ll2QQ7h , we determine the field-to- 
field variation in the blob number density. We then use 
a large volume cosmological N-body simulation to con- 
strain their host halo masses for the first time. 

In $2j we describe our narrow-band imaging survey and 
the data reduction procedures. $3] describes the selection 
of the Lya blob sample. In we present the Lya blob 
candidates (SJj, compare their s izes and luminosities 
to those of compact Lya sources ( §4.2)) , con strain their 
field-to-field variation and halo masses (S3]), and char- 
acterize their large-scale environment using the compact 
Lya emitters ( §4.4|) . In we discuss our conclusions. 
Throughout this paper, we adopt the cosmological pa- 
rameters Ho = 100/ikms -1 Mpc -1 , h = 0.7, = 0.3, 
and ft a — 0.7. All magnitudes are in the AB system 
(IOkdll974D . 

2. OBSERVATIONS AND DATA REDUCTION 

Using the MOSAIC-I and II CCD imagers on the 
KPNO Mayall and the CTIO Blanco 4m telescopes, we 
obtain deep narrowband images with a custom narrow- 
band filter (hereafter NB403 or NB). This narrowband 
filter has a central wavelength of A c « 4030A, designed 
for selecting Lya-emitting sources at z ~ 2.3. Its band- 
width of AAfwhm ~ 45A provides a line-of-sight depth 
of Az ~ 0.037, corresponding to 46.8 Mpc at z = 2.3 in 
the comoving frame. The CDFS, CDFN, and COSMOS 
survey fields have extensive ancillary data sets, includ- 
ing the deepest X-ray images for robust identification of 
AGN. 

We conducted the narrowband imaging observations 
over nine photometric nights between January 2007 and 
February 2009. The MOSAIC I and II cameras have 
eight 2k x 4k CCDs with a pixel scale of Of! 27 pixel -1 , 
leading to a sky coverage of 37' x 37'. We obtain deep 
narrowband images for four different pointings: CDFS, 
CDFN, and two COSMOS subfields (hereafter COS- 
MOS1 and COSMOS2). The total exposure time ranges 
from 7.2 to 10 hr, which consists of individual 20 or 30 
minute exposures with a standard dither pattern to fill in 
the gaps between the eight chips. The seeing ranges from 
l'/O to 1//3 depending on the fields. Table [1] summarizes 
our narrowband observations, including the central coor- 
dinates of the survey fields, total exposure times, survey 
areas, seeing, and survey depths. 
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Fig. 1. — Filter response profiles for the IVB403 narrowband (dot- 
dashed line) and broadband (17 and B) niters (solid lines) used in 
this study. The profiles are normalized to a maximum throughput 
of 1 and include the transmission of the atmosphere, telescope, 
camera optics, filter, and detector. ( Top) E SQ 2.2m 17 and B band 
filters from MUSYC (jGawiser et al. 2006a : b). (Middle ) KPNO U 
and Subaru B filters from the Hawaii-HDF survey (Capak et al. 
2004). (Bottom) C FHT u* and Subaru B filters from COSMOS 
(ICapak et al.ll20Q7h . 

To identify line emission objects requires that we sub- 
tract the continuum emission underlying the NB403 
bandpass. We estimate the continuum using existing, 
deep broadband (17 and B) images. For the CDFS, we 
use optical images from the Extended Chandra Deep 
Field- South dataset f rom the Mult i wavelen gth Survey by 
Yale-Chile (MUSYC; iGawiser et aTll2QQ6al MlFl. For the 
CDFN, we use the imaging produ cts from the Hawa ii 
Hubble Deep Field North survey (jCapak et al.l 12004ft . 
For COSMOS1 and COSMOS 2, we use 25 image tiles 
for each subfielcQ produced bv ICapak et al.l ((2007). For 
CDFS and CDFN, the broadband images have smaller 
sky coverage than our narrowband images, limiting our 
final survey areas (Tabled]). Table [2] lists details of the 
broadband images, including the filter properties, survey 
depth, seeing, and the instruments used. We also show 
the narrow and broadband filter transmission curves for 
all fields in Figure [TJ 

We reduce the narrowband image s with the IR AF ms- 
cred mosaic data reduction package (fVa ldes 1998) follow- 
ing th e procedures of the N OAO Deep Wide Field Survey 
team ((Jannuzi fc Devl[T999f ). The data are corrected for 
crosstalk between amplifiers and bias-subtracted. For 
flat-fielding, we use dome flats together with night-sky 

6 The original data were taken with the ESO MPG 2.2 m 
and Wide Field Imager ( WFI) by the ESQ Deep Public Sur- 
vey and COMBO - 17 teams ([A rnouts et al. 2001; Wolf et~aTll200l 
Erben et al. 2005; Hildebrandt et al. 2006). In this paper, we use 
the data products delivered by the MUSYC team. 

8 The tile numbers for the lower-left and upper-right corners of 
COSMOS1 and COSMOS2 are (17, 69) and (65, 117), respectively. 
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flats, which are median- combined from unregistered ob- 
ject frames each night. Satellite trails, CCD edges, bad 
pixels, and saturated pixels are masked. Th e astrometry 
is calibrated with the USNO-B1.0 catalog (jMonet et all 
2003) using the IRAF ccmap task. The individual im- 
ages are transformed to have the same pixel scales and 
world coordinate systems (WCS) as the reference broad- 
band images. For the COSMOS broadband images with 
finer pixel scale (0'/15 pixel -1 ), we resample them with 
the coarse MOSAIC-II pixel scale (Of! 27 pixel -1 ) to make 
the reference images. Finally, the projected images are 
scaled using common stars in each frame and stacked 
to remove cosmic rays. For flux calibration, we employ 
the 3-5 spectrophotometric standard stars observed each 
night to derive extinction coefficients and zero points for 
the NB403 magnitudes. Typical uncertainties in the de- 
rived zero-points are 0.02-0.04 mag. These uncertainties 
are added to the photometry errors in quadrature. 

In addition to the standard reduction procedures de- 
scribed above, we pay special attention to the "crosstalk" 
that occurs between CCDs (or amplifiers) sharing the 
readout electronics. The net effect of crosstalk is that a 
very bright or saturated source on one chip (or amplifier) 
produces echos or ghost images on the paired chips with 
proportional intensities, typically 0.02-0.2%, up to sev- 
eral counts. One can determine these proportionalities, 
known as crosstalk coefficients, using the science image 
itself, and thus remove echo images. In most circum- 
stances, this crosstalk correction works well and does not 
affect the detection and photometry of bright sources. 
Furthermore, dithered exposures often average out the 
ghost images in the individual frames when stacked to 
make the final products. 

However, residual ghost images from imperfect cor- 
rection can be easily confused with low surface bright- 
ness objects. In particular, for CTIO MOSAIC-II, the 
crosstalk images of a source at pixel (x, y) on one chip 
appears at the same pixel coordinates on the other chips 
(victim chips). Therefore, ghost images always appear at 
the same location relative to real sources in the dithered 
exposures, and they will be co-added in the final com- 
bined images, leading to false detections mimicking dif- 
fuse extended emission. 

To weed out false detections, we repeat the entire re- 
duction procedure doubling the crosstalk coefficients so 
that the ghost images are over-corrected and appear as 
negative counts. If a blob candidate is indeed an arti- 
fact arising from imperfect correction, it will show up as 
a negative image in the final combined image. Figure [2] 
shows an example in which one artifact image appears 
as a negative mirror image, whereas a real object is not 
affected by this crosstalk over-correction. Note that the 
artifact is as bright as the real blob candidate, demon- 
strating that special care is required to reject false de- 
tections for extremely low-surface brightness objects like 
Lya blobs. 

3. SELECTION OF Lya BLOB CANDIDATES 

To find Lya blob candidates, we construct photomet- 
ric catalogs in the NB403 narrowba nd and two broad- 
band s (U and B) using SExtractor (jBertin fc Arnoutsl 
1996|). First, we make "detection" images (NB+U+B) 
by adding the NB403 and broadband images after scal- 
ing them according to their signal-to- noise ratios (S/N). 
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Fig. 2. — Distinguishing between a false detection arising from 
an imperfect crosstalk correction of the CTIO MOSAIC-II image 
and a true extended Lya source. (From left to right) Original 
JVB403 images, NB403 images over-corrected with crosstalk coeffi- 
cients that are two times too large, and B band images. The ticks 
represent 10" intervals. Note that the residual due to the imper- 
fect crosstalk correction arising from the standard reduction (top 
left panel) resembles an Lya source that is spatially more extended 
than its continuum counterpart (top right panel). This artifact ap- 
pears as a negative mirror image in the over-corrected image (top 
middle panel), while the real object is not affected by the crosstalk 
over-correction (bottom middle panel). By comparing the original 
and over-corrected images, we are able to reject false detections, 
which is critical given the small number statistics of the blob sur- 
vey. 



After identifying sources in the "detection" images that 
have least 2 pixels that are 1.5a above the local sky, we 
run SExtractor in double-image mode on the JVB403, (7, 
and B images using these detection images. In other 
words, we first find the sources in the "detection" im- 
ages and then obtain photometry at their position in 
the JVB403, (7, and B images to construct three separate 
(one narrow- and two broadband) catalogs. We adopt 
Kron-like elliptical aperture magnitudes (i.e., MAG_AUT0 
in SExtractor) to derive photometric properties. Our use 
of the "detection" images ensures that 1) all the sources 
detected in either NB or the broadbands are included in 
our catalog and 2) the elliptical apertures of the more ex- 
tended sources in the (NB+U+B) image are large enough 
to include all the light from both the NB and broadband 
images. Note that this choice of photometric aperture 
is different from the classical Lya emitter searches that 
adopt a small circular aperture (a few x FWHM of see- 
ing) to detect fainter sources. 

The selection of Lya blob candidates from the NB and 
broadband photometric catalogs consists of two steps: 1) 
selection for line (hopefully, Lya) emitting objects with 
large line equivalent widths and 2) selection for spatially 
extended objects with a larger angular extent in line 
emission than in the broadbands. In other words, we 
define a "blob" as an object whose Lya emission above 
a certain surface brightness threshold is more extended 
than its stellar continuum, thus representing light from 
the intergalactic medium. 

First, we choose candidates by requiring that they are 
detected above the completeness limits of the JVB403 
images (JVB403 < 24.5 mag). All candidates must 
have observed- frame equivalent widths larger than 100 A 
(EW r est > 30.3 A), corresponding to (UB - NB) > 
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Fig. 3. Color - magnitude (UB-NB403) versus JVB403 dia- 
gram for all the sources detected in our four survey fields: CDFS, 
CDFN, and two COSMOS subfields. Here UB represents the AB 
magnitude of the average continuum flux density within the JVB403 
band estimated from the U and B band images. Right and top axes 
show the corresponding equivalent widths in the observed frame 
and the NB fluxes, respectively. We select line-emission objects 
with the criteria NB < 24.5 (solid line) and (UB-NB403) > 1.2 
(dashed line; EW Q b s > 100A). Open circles represent the final Lya 
blob candidates obtained from Fig. [4] There are more open circles 
than final blob candidates in the CDFS panel because a few blob 
candidates include multiple sources. 



1.2, where UB represents the AB magnitude of the av- 
erage continuum flux density within the NB403 band 
(fcont) estimated from the U and B band images, UB 
= -2.5 log(/^ t ) -48.60. We estimate f™ t and the line 
flux (Fu ne ) of these objects using the following relations: 
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where Fjj, Fb and F^b are the total flux in each filter 
derived from the 17, B, and NB403 magnitudes, respec- 
tively. A Ac/, AA#, and AXnb represent the band- widths 
of the U", B and NB filters, respectively. f^ ont and fcont 
represent the average flux density of galaxy continuum 
within the U and B bands, respectively. €b and ejj are 
the correction factors that we use to remove the NB light 
from each broadband when estimating fj/ ont and fcont- 
Figure [3] shows the (UB — NB) color as a function of NB 
magnitude for all objects detected in either the narrow 
or broad bands within our survey area. 

In CDFS, COSMOS1, and COSMOS2, bright (18 < 
NB < 22) sources, mostly stars and nearby galaxies, have 
an average (UB — NB) color that is consistent with zero 
(-0.llzb0.18, 0.13±0.17, and 0.12±0.12, respectively). 



In the CDFN, however, they have an average color of 
(UB — NB) ~ —0.24 ± 0.14 mag, which suggests absorp- 
tion features in our narrowband that are very unlikely. 
Because this —0.24 mag offset disappears when we adopt 
continuum measurements (in u and g) from the Sloan 
Digital Sky Survey, we apply a small c orrection (0.2 mag) 
to the CDFN broadband magnitudes (jCapak et al.l l2004) 
such that (UB—NB) ~ for the bright sources. Although 
we cannot explain the appa rent problem with t he CDFN 
broadband photometry in iCapak et al.l ([2004), the se- 
lection of our blob sample is unaffected by the applied 
correction. 

At our survey redshift, the only possible interlopers are 
nearby [O II] A3727 emitters at z ~ 0.08. However, such 
objects rarely h ave equivalent wid ths larger than 100A in 
the rest frame (|Hogg et al.lll998f ). Therefore, we expect 
that the contamination of our z = 2.3 Lya source catalog 
by nearby star forming galaxies is minimal. 

Second, we identify those line-emission selected objects 
that are more spatially extended in Lya than their con- 
tinuum counterparts (Fig. 2]). This selection d efinit ion is 
the same as that adopted by Matsuda et al.l (12004) and 
somewhat different than that of iSaito et al.l (|2006), who 
select spatially extended objects by requiring the FWHM 
in their intermediate-band to be larger th an that in the 
broad band or continuum image (see also Nils son et al.l 
[20Q9h . 

We measure the spatial extent of the Lya emission in 
the continuum-subtracted images. After registering the 
NB403 and broadband images ( U and B) at the sub-pixel 
level and matching their seeing, we construct continuum- 
subtracted NB403 images by applying the relations in 
Eq. (pQ) in 2-D. We measure the isophotal area of the 
emission region by running SExtractor with a thresh- 
old of 5.5 x 10 -18 ergss -1 cm -2 arcsec -2 , which corre- 
sponds to 3cr, 1.8a, 2.2a, and 2.1a above the local sky in 
the CDFS, CDFN, COSMOS1, and COSMOS2 fields, re- 
spectively. This measu rement threshold is ~ 2.5 x higher 
than that adopted by iMatsuda et all (|2004j ). However, 
because our survey redshift (z = 2.3) is lower than theirs 
(z = 3.1), we gain a factor of ~ 2.4 in surface brightness, 
thus achieving equivalent surface brightness (8.6 x 10 5 Lq 
kpc -2 compared to their 8.2 xlO 5 L kpc -2 ). We also 
estimate local sky background using a fairly large back- 
ground mesh size (~ 60"x60"), so as not to mistak- 
enly subtract the extended Lya emission as a local back- 
ground. 

Measuring the size of a low surface brightness feature 
is always subject to the noise and filtering. After testing 
the various smoothing filters in SExtractor, we adopt 
a 5x5 pixel 2 Gaussian kernel with FWHM = 2 pixels 
to compromise between signal-to- noise (S/N) and over- 
smoothing. We choose this truncated smoothing kernel 
such that it can enhance S/N in the low surface bright- 
ness wings while not spreading the bright core into the 
outer part. We note that the measured isophotal area de- 
pends strongly on the choice of the filters and recommend 
specifying the smoothing kernel and detection threshold 
when reporting blob size to allow a direct comparison be- 
tween different samples. A ppendix 1X1 shows wha t sizes 
and fluxes the 35 blobs from lMatsuda et all (|2004t ) would 
have if those blobs were m easured followinjg; our proce- 
dures. In our survey, the IMatsuda et al.l (|2004f ) blobs 
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Fig. 4. — Distribution of isophotal areas (Ai so ) and Lya lumi- 
nosities of Lyai-emitting sources (red circles) in the CDFS, CDFN, 
and two COSMOS subfields. Thick solid lines represent the selec- 
tion criteria for the final Lya blob candidates (filled circles): A- lso 
> 10 D ,f and distinguishable from point sources, i.e., deviating from 
the simulated point source locus of small grey dots by > 4cr. Simu- 
lated point sources mix with extended Lya-emitting sources below 
our selection limits in A[ so and L-^ ya , making it difficult to dis- 
tinguish blobs from point sources there. The horizontal dashed 
line represents the division (A- lso = 16 □") between the brightest, 
largest blobs and the other blob candidates. The former are dis- 
covered only in the CDFS. The green squares in the CDFN, COS- 
MOS1, and COSMOS2 panels represent the A[ so and L^ ya that 
the 16 CDFS blobs would have if they we re ob served at the seeing 
and depth of each of the other fields (see M4.3)) . The six brightest, 
largest blobs (L Lya > 1.5xl0 43 ergss -1 and A iso > 16 in the 
CDFS should have been detected in the other fields. Therefore, 
we conclude that the observed strong fiel d-to- field variation is real 
and not due to observational biases (see 34.30 . 



would be smaller, and the smallest blobs would have an 
isophotal area of A{ so ~ 10 Therefore, we adopt this 
value as the size lower-bound when selecting our blob 
sample. 

One of the potential problems in detecting a blob is 
contamination by point sources. To quantify this effect, 
we place artificial point-sources with a range of luminosi- 
ties (i^Lya = 10 42 - 10 45 ergss -1 ) into the sky regions 
and measure their sizes and fluxes in the same manner as 
for the extended sources. We then determine the isopho- 
tal area versus luminosity limits above which extended 
and point sources can be differentiated. 

Figure |4] shows the distribution of the angular sizes 
and line luminosities of the Lya-emitting blob candidates 
assuming that they are all at z = 2.3. The open and filled 
circles represent the line-emitting objects selected using 
our line-emission criteria (Fig. [3]), and the gray dots show 
the relation between size and brightness for the artificial 
point sources. To choose the final blob candidates, we 
select objects with isophotal areas larger than 10 U" that 
lie more than 4a above the L\, ya -Ai SO relation defined by 
the point sources. Below these limits, extended and point 
sources mix, and the sizes of blobs cannot be measured 
reliably, as explained below. 
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Fig. 5. — Cut-out images of the 16 Lya blobs that we dis- 
cover in the Extended CDFS. Images from left to right: U, 
IVB403, continuum-subtracted Lya line (A c — 4030A), B, and 
R band. The ticks are spaced in 10" intervals. The overlayed 
(yellow) contours represent a Lya: surface brightness of 4xl0 -18 
ergss -1 cm -2 arcsec -2 , our 2.2a detection limit in the narrowband 
images. In each case, the Lya emission is more extended than the 
broadband counterpart. All five candidates that we have now ob- 
served spectroscopically are confirmed as blobs at z = 2.3 (Yang 
et al., in prep.). 

Because the chosen isophotal threshold is comparable 
to the rms sky noise (1.8a - 3a depending on the field), 
we test how reliably we can measure the spatial extent 
of the blob candidates. We cut a small (101x101 pixel 2 ) 
section around each candidate from the line-only image, 
filter it with a smoothing kernel^ place each postage 
stamp into ~ 1000 empty sky regions, and extract the 
sources with SExtractor in the same way as for the real 
data. Then we check how often and accurately the blob 
sizes are recovered from these simulated images. The 
recovery fraction (/ recv ) represents how frequently an ar- 
tificial blob is recovered with a size larger than 10 
The size error becomes comparable to the measured size 
below our sample selection limits, and f reC v drops from 
>90% to ^50% at the selection boundary for the bright- 
est blob candidates. 

Our "blob" definition requires their Lya emission to 
be more spatially extended than their stellar continuum 
emission. Thus, their extended light represents the in- 
tergalactic medium instead of the galaxies. Our choice 
of size cut is not due to any discontinuity, and is lim- 
ited only by the ground-based seeing. Therefore, higher 
spatial resolution images would likely detect fainter or 
more compact blobs. Stacked images of LBGs (Ly- 
man Break Galaxies) d o reveal faint extended emission 
(|Havashino et al.ll2004D . 

4. RESULTS 
4.1. Discovery of Lya Blobs 
The total area of our survey is ~ 1.2 deg 2 . Our Lya 
blob (LAB) selection criteria, EW Q bs > 100 A and A iso 
> 10 □" above the surface brightness threshold of 5.5 x 



10 



-18 



ergss cm 2 arcsec , yield a total of 16, 5, 3, 



1 We use a 5x5 pixel 2 convolution mask with FWHM = 2 pixel. 
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Fig. 5. — Continued. 



Fig. 5. — Continued. 




Fig. 5. — Continued. 



and 1 blob candidates in the CDFS, CDFN, COSMOS1, 
and COSMOS2 fields, respectively. We list their basic 
properties in Table [3) coordinates, Lya luminosity, and 
isophotal area. 

Figures [SHU show postage-stamp images of all 25 blob 
candidates in the ATB403 band, continuum-subtracted 
Lya line emission, and three broadbands (UBR) over- 
layed with the Lya contour corresponding to a surface 
brightness of 4 x 10 -18 ergss -1 cm -2 arcsec -2 . 

Although we select blob candidates quantitatively, 
follow-up visual inspection indicates that the spatial 
extents of four blobs are not clearly larger than the 
local PSF: CDFN-LAB02, CDFN-LAB03, COSMOS- 
LAB01, and COSMOS-LAB02. Therefore, we flag them 
as "marginal". Excluding them does not affect our 
conclusion regarding the large field-to-field variation in 
the blob number density and, in fact, makes that case 
stronger. The largest blobs (more than 16 □"; the dashed 
line in Fig. [4]) are the most robustly extended and are 



comparable to those identified in other surveys (e.g., 
iMatsuda et al.l 12004} ). Therefore, we treat them sepa- 
rately as the "bright/large" subset of the entire blob 
sample in the clustering analysis in §4.31 

The blob candidates have a wide range of sizes and 
line luminosities: A{ so = lOD"- 60 U" and L\, yot = 
0.7-8 x 10 43 ergss -1 . They show diverse morphologies 
ranging from compact (the south-west clump of CDFS- 
LAB01) to diffuse (CDFS-LAB05) to highly elongated 
(CDFS-LAB08 and 16). We do not find any bubble-like 
structures that might be associated with superwinds like 
those in a few of the blobs identified by Mats uda et al.l 
([2001 . 

We are following up the entire blob sample spectro- 
scopically. To date, we have observed five blobs (CDFS- 
LAB01, 02, 04, 10, 14) of the 16 in the CDFS. We confirm 
all five spectroscopically with their Lya and/or Ha lines. 
The details of our spectroscopic campaigns are presented 
in a forthcoming paper (Yang et al., in prep.). 

The properties of continuum objects associated with 
blobs provide valuable clues to the source of blob emis- 
sion (Yang et al., in prep.). For e xample, the bl o b dis- 
covered in the CDFS at z = 3.1 bv lNilsson et al.l (|2006f ) 
is not associ ated with any con tinuum source within ~3" 
(but see also lGeach et al.ll2QQ9h . Thus, it is possible that 
this blob is powered by cooling radiation or cold mode 
accretion. Although most blob candidates in our sample 
have clear continuum source counterparts in the rest- 
frame UV images, two blobs (CDFS-LAB04 and CDFS- 
LAB05) do not. We are, however, able to identify con- 
tinuum sources within these blobs in deep, rest-frame 
UVHST images (GEM S and GOODS-S: lRix et al.ll200l 
iGiavalisco et al.l l2004f ). Preliminary inspection of the 
multiwavelength images also reveals bright IR or X-ray 
sources in these blobs. Interestingly, we often find multi- 
ple continuum sources in a blob, notably in LAB02 and 
LAB03. Our initial examination of the HST images of 
the 16 Lya blobs in the CDFS often resolves these sources 
as galaxies, suggesting that star formation and/or nu- 
clear activity might play a role i n producing the Lya 
emission (e.g.. iColbert et al. 2006). This apparent clus- 
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Fig. 6. — Same as Figure [5] but for the five Lya blob candidates 
in the CDFN field. 



tering of sources within some blobs may indicate that 
blobs are the progenitors of groups or clusters of galax ies 
today. We test this possibility quantitatively in §4.31 

4.2. Continuous LAE-to-LAB Sequence 

There is no evidence of a discontinuity between the 
properties of our blobs and unresolved, compact Lya 
emitting (LAE) galaxies (i.e., the sources below the 
solid lines in Fig. |4]). The smoothness of the LAE- 
to-LAB transi t ion, w hich has also been observed by 
iMatsuda et al.l ([2004), suggests that these two types of 
Lya sources may not be distinct and that whatever mech- 
anism or mechanisms power blobs work over a wide range 
of luminosity and spatial extent. Understanding the ori- 
gin of the extended Lya emission requires us to probe 
the host galaxy properties — including the stellar mass, 
star formation rate, and size of the star- forming region 
— along the emitter-to-blob sequence. If extended star 
formation or AGN are responsible for powering the emis- 
sion, we would expect that the properties of blobs (A[ so , 
Lhya) are correlated with either star formation rate or 
X-ray luminosity of the host galaxies. For example, 
IGeach et all (2005) argue that a correlation between the 
Lya and bolometric (actually FIR) luminosity (although 
weak) suggests that the interaction of an ambient halo 
of gas with a galactic-scale superwind is responsible for 
the majority of LABs. We discuss the multi-wavelength 
properties of galaxies within or near the blobs, consider- 
ing them as possible energy sources, in a separate paper 
(Yang et al., in prep.). 

4.3. Significant Differences in Blob Counts per Field 

Surprisingly, most blobs (16/25) and all eight of the 
brightest, largest blobs (A iso > 16D", L Lya > 1.5xl0 43 
ergss -1 ) lie in only one of the survey fields: CDFS. Be- 
cause the depth and seeing of the CDFS images are also 
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Fig. 7. — Same as Figure(5] but for the three Lya blob candidates 
in the COSMOS 1 subfield (top rows) and the one candidate in the 
COSMOS2 subfield (bottom row). 



superior to the other fields, we first need to verify that 
this field-to-field variation of the blob number density 
does not arise from selection effects. For example, al- 
though the separation between the fainter /smaller blob 
candidates and the simulated point-sources is distinct in 
the CDFS, the separations in other three fields are less 
clear. 

To confirm that the field-to-field variation of the blob 
population is not due to the different seeing and sur- 
vey depth, we simulate how the blob candidates in the 
CDFS field would look if observed in the other three 
fields. Using the continuum-subtracted postage stamp 
(27"x27") images of the CDFS blob candidates, we first 
degrade their seeing to that of the CDFN, COSMOS1, 
and COSMOS2 fields by convolving kernels derived from 
PSF images. We rebin the images (if pixel scale is differ- 
ent), add Poisson noise, and place them into empty sky 
regions in the continuum-subtracted narrowband images 
of the other fields. We then measure the blob sizes and 
luminosities in the same way as described in ^3] We re- 
peat this experiment ~ 1000 times to derive the range of 
recovered luminosities and sizes. 

We show the isophotal areas and Lya luminosities 
(Ai SO ~ ^Lya) from this recovery test as green squares 
in the CDFN, COSMOS1, and COSMOS2 panels in Fig- 
ure |U If we estimate the number of possible detections 
as N = Xl/ recv >50% /recv? the number of CDFS-like blobs 
that should have been detected as bright/large blobs with 
A iso > 16 U" and L Lya > 1.5 xlO 43 ergss -1 is 6.9, 6.0, 
and 6.0, respectively, for each of the three other fields. 
However, we do not find such blobs in other fields, so we 
conclude that the observed field-to-field variation is real 
and that the minimum variation in the number counts 
of the largest blobs (A{ so > 16D") is at least iVbiob = 

(TVcDFS, N CDFN i ^COSMOS!, ^COSMOS2) = (6, 0, 0, 0). 

Throughout the rest of this paper, we refer to this as the 
"bright/large blob sample" and use it as the default for 
analyzing blob statistics. When we lower the criteria for 
the blobs to A iso > 10D", we expect 14.0, 13.9, and 13.7 
blobs in the other three fields, respectively, whereas we 
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find 5, 3, and 1 blob candidates. In this case, the ob- 
served contrast between the four survey fields is iVbiob = 
(14, 5, 3, 1). We subsequently refer to this as the "entire 
blob sample". We consider the bright/large sample and 
entire sample field-to-field variations in Abiob? which are 
corrected for the different survey conditions, in deriving 
the average number density and its variance in the fol- 
lowing paragraphs. The blob statistics are summarized 
in Table gj 

Could this observed field-to-field variation arise solely 
from the statistical uncertainty? Before we proceed to a 
detailed analysis, we have to rule out the possibility that 
this result arises solely from Poisson (shot) noise. For the 
hypothesis that the surface density of the bright/large 
Lya blob sample is 2-20 blobs deg -2 or ~ 0.5 - 5 blobs 
per survey field, we calculate the probability of observ- 
ing non-detections in three survey fields and detections 
more than six blobs in any of the fields. We are able to 
rule out a uniform distribution with at least 99.97% con- 
fidence (~3.6<j). For the entire blob sample, a uniform 
distribution is excluded at the 3.7a level. In §4.3.3| we 
compare the observed field-to-field variations with those 
derived from Poisson statistics and cosmological N-body 
simulations, respectively. We show that the latter better 
reproduces the observations. 

4.3.1. Quantifying Blob Field-to- Field Variance 

In this section, we estimate the field-to-field variation 
in the number density of Lya blobs from the observed 
number statistics, iVbiob- According to ACDM cosmol- 
ogy, the number density and variance of a galaxy popula- 
tion are not entirely independent properties, but a func- 
tion of halo mass. Here we treat them as independent 
parameters and aim to measure them as observables by 
adopting a simple analytic approximation of the under- 
lying fluctuations in blob number density arising from 
large-scale structure. This method has the advantage 
that the blob number density and variance can be derived 
over a wider range of parameter space than sampled by 
simulations. We then compare these properties with the 
predictions from the simulations to obtain constraints on 
the halo mass ( §4.3.3p . As will be shown in the following 
section ( §4.3.2|) . one must consider the field- to-field vari- 
ations to correctly estimate the uncertainties in number 
density and to thus compare blob statistics across differ- 
ent surveys. 

The simplest way of quantifying the field-to-fi eld vari- 
ation, a 2 , is to adopt the relation (Peebles 1980, §36): 

2 (A" 2 ) - (AQ 2 _ J_ 

(N) 2 (N) ' K ) 

where N is the number of blobs per 0.25 deg 2 , the typical 
area of each of our four fields, and a 2 is the fractional 
variance corrected for Poisson noise. The number den- 
sity n can be derived directly from A^ by dividing it by 
the survey volume: n = 8.6 x 10 -6 N Mpc -3 . Hereafter, 
we use surface density to mean surface number density. 
Thus (j 2 represents the fractional uncertainty in the ob- 
servational estimate due to finite survey volume. For the 
bright/large sample, we obtain (N) = 1.5 and a v ~ 1.5 
(150%). For the entire blob sample, (N) = 5.75 and a v 
~ 0.76, or 76%. Because of our small number statistics, 
we choose to adopt a more sophisticated method to bet- 



ter understand the possible range of blob number density 
and its variance. 

To quantify the field-to-field variance of the blob pop- 
ulation, we calculate the posterior probability for a 2 and 
an average surface density per 0.25 deg 2 (N) given our 
observation (D) of (6, 0, 0, 0) blobs in the bright/large 
sample: 

p(a v ,N\D)(xprob(D\a v ,N) p(a v ,N). (3) 

First, we assume that the field-to-field variance follows 
the log- normal distribution: (N/N) ~ Log-N(0, <Jl N ); 
in other words, log(N/N) follows a normal distribu- 
tion A'(0,cr 2 jN ). Here, <Tl N is the variance of the log- 
normal distribution and is related to the actual variance 
by a 2 = exp(aL N ) — 1. Unlike a Gaussian distribution, 
the log-normal distribution does not allow negative val- 
ues for N and naturally introduces a skewness into the 
distribution. When ^ < 1, the log- normal distribu- 
tion is similar to Gaussian, but the distribution becomes 
skewed toward zero as a v increases, effectively mimick- 
i ng the dark matter fluctuations at the high mas s end 
(jColes fc Jones! [1991 iBernardeau fc Kofmanl fl995l ) . We 
choose the log- normal distribution for the simplicity here, 
but any reasonable functional form capable of represent- 
ing this skewness can be used. 

Second, for a given set of (a v ,N), we calculate the 
probability, prob(D\a Vl AT), of finding six blobs in one 
field and none in three other fields assuming that the 
observed number of blobs follows Poisson statistics with 
a mean of N. We adopt logarithmic priors for both a_ v 
and N, which indicates p(cr v ) oc l/cr v or p(N) oc 1/N, 
implying that the scale of N and a v is unknown, i.e., that 
the priors are uniform in logarithmic bins. We consider 
a range of 0.1 < A" < 30 and 0.1 < a v < 10, i.e., 10% to 
1000% field-to-field variance. We also test other priors 
including (1) a linear prior for both a v and A" and (2) a 
logarithmic prior for a v and linear prior for A", but the 
choice of prior does not affect our conclusions. 

Figure [8] shows the posterior probability distribution 
of the average surface density and variance for the 
bright/large blob sample (left) and the entire sample 
(right). For the bright/large blobs, the posterior favors 
high variance (a v ^> 1) as expected. The confidence 
regions are not closed, allowing us to put only a lower 
bound on a v . The lower limits are a v > 1.45 and a v > 
0.57 for the la and 2a confidence levels, respectively, for 
the joint distribution (i.e., we attempt to constrain both 
N and a v at the same time). The la, 2a and 3a confi- 
dence regions are determined such that the correspond- 
ing contours enclose 68.3%, 95.5%, and 99.7% of total 
posterior probability. The 3a lower limits are not con- 
strained because the priors become too high as a v — > 0, 
thus the posterior probability is determined more by the 
input priors rather than by the data themselves. 

We also calculate the confidence interval for each pa- 
rameter (a v or AT") by marginalizing out the other pa- 
rameter in the posterior probability distribution: e.g., 
p(N\D) = J p(a Vl N\D)da v . This interval then repre- 
sents the estimate of one parameter independent of the 
other. In £|4.3.2[ we compare the marginalized surface 
density estimates, regardless of the variance, with those 
obtained from past observational studies. In Table [5j 
we show the blob surface density (N) and variance esti- 
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Fig. 8. — Posterior probability distribution of the fractional variance (a v ) and the average surface density per field (N) for the bright, 
large blobs (left) and the entire blob sample (right). We adopt logarithmic priors for both a v and N. The la, 2a and 3a confidence regions 
are determined such that the corresponding contours enclose 68.3%, 95.5%, and 99.7% of the total posterior probability. The bright, large 
blobs have fractional variances of at least 1.5 (150%) at la and 0.6 (60%) at 2a regardless of the estimated number density, confirming 
that blobs are highly clustered. For the entire blob survey, the most probable estimate is (N, a v ) ~ (4.8, 0.9). 



mates obtained from the above. The field-to-field vari- 
ation of the bright/large blobs is stronger than 150% 
and 60% at the la and 2a confidence levels, respectively. 
For the entire blob sample, the most probable parame- 
ters are (TV, a v ) ~ (4.8, 0.9) with large uncertainties of 
0.26 < a v < 5.2 at the la level. The large variance im- 
plies a highly skewed distribution and thus reproduces 
the observed strong field- to-field variations. 

Observationally, a v ~ 1.5 (150%) is much larger than 
for compact Lya emitters (LAEs) at z = 3-5 obtained 
from narrowband imaging over volumes comparable to 
our survey. For example, we estimate a LAE variance 
of ~ 20-30% at z^3 .1 from the five 0.2 deg 2 s ubfields 
of lOuchi et all (|2008f ). From lShiova et al.l (|2009D . we es- 
timate that the LAE variance at z ~ 4.86 is likewise 
~ 30% when their contiguous survey area (~ 2 deg 2 ) 
is divided into 0.25 deg 2 subfields (their Table 2). This 
~30% variance is enough to produce the factor of 2 differ- 
ence among the LAE number densities observed in their 
fields. Note that the ~30% LAE variance is uncorrected 
for Poisson noise, which would lower it. Therefore, the 
even larger variance of blobs suggests that they lie in 
halos more massive than those of LAEs. We put the 
first constraints on the dark matter halo mass of blobs 
in §4.3.31 by comparing our blob statistics with N-body 
simulations. 

The strong variation in blob counts from one 0.25 deg 2 
field to another is consistent w ith the discovery of a close 
pair of blobs by lYang et aLl (j2009[ ). Those blobs were 
among only four detected in our shallower, but larger 
(4.8 deg 2 ) survey of the NOAO Bootes field. They are 
separated by only ~70". If that survey had been con- 
ducted by mosaicing 30'x30' fields like those sampled 
by our imager here, most fields would not contain the 
pair. On a related n ote, the clustering of 35 blobs in the 
SSA22 overdensity (Mat suda et al.ll2004f ) was identified 
during a follow-up of two giant blobs originally found by 



iSteidel et all (|2QQQh . 

4.3.2. Blob Number Density and Luminosity Function 

The strong field-to-field variation of the Lya blobs 
presents challenges for the measurement of their number 
density and luminosity function (LF). In this section, we 
compare our LFs derived from each survey field with each 
other and with previous studies. We demonstrate that a 
large volume survey and/or multiple pointings are crit- 
ical to constrain the blob LF. For the rest of paper, we 
adopt the marginalized number densities from the pre- 
vious section, which give us the number densities of n 
= l.OtojxlO- 5 Mpc- 3 (from N = 1.2 per 0.25 deg 2 ) 
for the bright/large blob sample and n = 4.ll^gXl0 -5 
Mpc -3 (from N = 4.8 per 0.25 deg 2 ) for the entire blob 
sample. 

First, we compare the nu mber density of bright/large 
blobs derived from §4.3.11 with previous measurements 
(Figure [9j). To make a fair comparison with other 
work, we need to cut the bright /lar ge sample t o sat - 
isfy the selection criteria used by Yan g et all ([2009) 
when they previously examined the blob number den- 
sities among different samples spanning z = 0.8 to 6.6. 
Those criteria were: L\, ya > 1.5 xlO 43 ergss -1 and A lso 
> 25 U" above a surface brightness threshold of 5xl0 -18 
ergss -1 cm" 2 arcsec -2 under ~1.8" seeing. We refer 
readers to I Yang et"aTI (j2009[ ) for the details about the 
estimates at each redshift. 

Only four blobs (CDFS-LAB01, 02, 03, and 04) 
from the s a mple of six bright/large blobs satisfy the 
lYang et"aTI (j2009[ ) criteria, so we scale the bright/large 
counts down by a factor of 4/6, obtaining a surface den- 
sity of 0.77±J-j per 0.25 deg 2 fields or a number density 
of n = 0.66lo 2 xlO -5 Mpc -3 (large, open square in 
Fig. [9]). The error bar takes into account the field-to- 
field variation as well as the Poisson noise. This average 
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Fig. 9. — Number density of Lya blobs at different redshifts. 
The filled and open squares are from our two narrow-band imag- 
ing sur veys at z = 2.3: the ultra-wide Bok+90Prime survey in 
Bootes ( Yang et al. 2009) and the narrower, but deeper, NOAO- 
4m survey (this work), respectively. Th e filled dots are the numbe r 
density estimates from Matsuda et al. (2004); Saito et al. (2006); 
IQuchi et all (|2009h ; Keel et al. ( 20®. All the number densities 
plotted here are scaled to match the shallowest survey (Bootes 
field) with ~1.8 7/ seeing and an rms sky background of 5xl0 -18 
ergs s _1 cm -2 arcsec -2 . The upper or low er limits on several of the 
points from the literature are explained in Yang et al. (2009). The 
four small squares represent the individual number density esti- 
mates for each of our four survey fields, demonstrating that a large 
volume survey is required to overcome the strong field-to-field vari- 
ation of the blobs. Note that these data points are shown slightly 
shifted for clarity, but they are all at z = 2.3. 



number density at z = 2.3 derived from the f our NOAO- 
4m su rvey fields is thus consistent with the lYang et"a"T1 
(2009) measurement at the same redshift (n = 0.25=b0.10 
xlO -5 Mpc -3 ; large, filled square) obtained from our ~ 
4x larger volume survey of Bootes. 

To illustrate the uncertainties arising from the field- 
to-field variance, we show the bright/large number den- 
sity for each of the four fields individually. The small 
squares represent n = 3.2 (±1.6) xlO -5 Mpc -3 for the 
CDFS, and three upper limits, n < 0.90, 0.62, 0.60 
xlO" 5 Mpc" 3 for CDFN, COSMOS1, and COSMOS2, 
respectively. The number den sity in the CDFS field is 
consistent with that found by iMatsuda et al.l (|2004f ) in 
the SSA22 field (5.8±2.1 x 10~ 5 ) within the uncertain- 
ties, suggesting that blobs at z = 2 .3 in the CDFS also 
lie in an overdense r egion (see also iPalunas et al.l 120041 : 
iPrescott et al.l [2QQ8h and occupy high mass halos. On 
the other hand, the blob number densities for the other 
three fields are lower than for the CDFS and the SSA22 
field, demonstrating that characterizing the evolution in 
n with z requires surveys large enough to overcome the 
field- to-field variance. Because such measurements are 
not yet available at other redshifts, it is not possible to 
constrain n(z) at this time. 

Note that we are not able to apply the 
(2009) sele ction criteria to the higher- z 
(2006§ and IQuchi et al.l (j2009f ) blob samples because 



those authors use se l ection methods different than ours 
and IMatsuda et al.l (|2004f ). Therefore, we plot these 
higher-z estimates as upper and lower limits, respec- 
tively. Unlike the comparison betwee n the z = 2 . 3 and 
z = 3.1 samples , comp arison of the iSaito et"aTI ([2006) 
and IQuchi et al.l (j2QQ9h points with our results is per- 
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Fig. 10. — Cumulative Lya luminosity function (LF) for 
Lya blobs. The open squares (red) and diamonds (blue) rep- 
resent the LFs derive d from the CDFS (this paper) and SSA22 
(Matsuda et al. 2004) fields, respectively, which each have a sur- 
vey area of ~0.25deg 2 . The filled dots represent the average LF 
measured from all four of our survey fields (CDFS, CDFN, COS- 
MOS1, and COSMOS2). Data points are slightly shifted relatively 
to each other for clarity. While the slope and normalization of the 
CDFS and SSA22 LFs roughly agree, the normalization of the av- 
erage LF is ~ four times lower, suggesting that one should account 
for the strong field-to-field variation of Lya blobs when calculating 
LFs. 



ilous. We plot all these points together only to summa- 
rize the state of blob surveys. 

In addition to the blob statistics obtained from 
narrow- or intermediate-band imaging techniques, we 
show the number density estimate (n < 0.5 xlO -6 
Mpc -3 ) from GALEX slitless spectroscopy of t wo galaxy 
(super) clusters at z - 0.82 ([Keel et al.l l2009D . This n, 
even though calculated from overdense regions, is well 
below that of the overdens e CDFS and [Matsuda et al.l 
(|2004f ) fields, supporting the lKeel et all (|2009f ) claim that 
Lya blobs might be high redshift phenomena. Additional 
surveys at z < 1 are required to confirm this result, as 
we do not know whether the variance in n at z = 2.3 
persists at lower z. 

Second , we compare our blo b luminosity function with 
that of IMatsuda etHI (j2004f ) in Figure [TUl The simi- 
lar depth and selection criteria of the two surveys makes 
the direct comparison meaningful (see Appendix [A]) . For 
each of our blobs, we use the total Lya luminosity (Ta- 
ble [3]), which our recovery test demonstrates is reliable. 
The isophotal flux could be biased due to the different 
signal-to-noise and filtering kernels used in the detection 
procedure. We construct an average LF as follows. In 
each Lya luminosity bin, we divide the number of blobs 
found in all four of our survey fields by the total sur- 
vey volume. We also plot the LF obtained from only 
the 16 blobs in the CDFS. We do not apply a complete- 
ness correction using the recovery fraction f recv in Ta- 
ble |3l because we also do not know the completeness of 
IMatsuda etHI (|2004f l sample. In Table El we list our 
blind s urvey LFs as well as t hat of the SSA22 overdense 
region (jMatsuda et al.ll2004fl . 

T he slope and normal ization of the CDFS (squares) 
and lMatsuda et aTl (j2004f ) (diamonds) LFs agree roughly. 
The apparent discrepancy in the faintest bin (log Z, LyCK = 
42.8) is likely due to incompleteness in our survey. We 
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also do not find any blobs as bright as those in their 
brightest bin (log L\, ya = 44.0). Not surprisingly, the nor- 
malization of our "average" luminosity function (filled 
circles) is ~ 4x low er than those of the CDFS and 
iMatsuda et al.l j2004) LFs. Once again, we see that com- 
parisons among surveys are difficult without knowledge 
of the blob clustering strength and that large volume 
surveys are required to overcome the strong field-to-field 
variance of the blob counts. The same conclusion is 
reached fr om the on-going Su baru survey of blobs (see 
Fig. 13 of iGoerdt et al.l 12009. Y. Matsuda in prepara- 
tion). 

The clustering of blobs provides a means to discrimi- 
nate among models for the origin of the extended Lya 
emission, including photoionization by AGN, outflows 
due to intense star formation, and cooling radiation. 
This analysis is beyond the scope of this paper, but we 
note that any blob-producing mechanism must repro- 
duce both the observed LF and its field-to-field varia- 
tion. For example, it would be interesting to ascertain 
whether AGN or sub-millimeter galaxies (SMGs) have 
clustering strengt hs sim ilar to blobs. Along these lines, 
iDijkstra fc Loebl (j2009f ) claim that Lya blobs are cool- 
ing radiation arising from cold streams falling onto the 
embedded galaxies and that a strong variance naturally 
arises fro m the underlying v ariation of dark matter halos 
(see also IGoerdt et al.l 2009). 

4.3.3. Blob Halo Masses 

Because our blob survey is blind and sufficiently large 
to determine the blob number density and its variance, 
we have a unique opportunity to constrain the properties 
of the dark halos in which the blobs reside, and thus to 
understand what these mysterious objects have evolved 
into today. For example, based on the blob number den- 
sity and the discovery of the blob pair i n the NOAO 
Deep- Wide Bootes field, lYang et all j2009) suggest that 
blobs are sites for the formation of the brightest galax- 
ies in rich galaxy clusters. However, the small number 
statistics of that study precluded constraining the mass 
of blob halos. Although our current survey statistics are 
still not large enough to directly measure the clustering 
of blobs via correlation function analysis, we can use n, 
a v , and a cosmological N-body simulation of ACDM uni- 
verse to identify the most likely halo mass occupied by 
blobs. 

We first select a dark halo (DM) mass in which blobs 
could reside. All halos above this minimum mass M m i n 
have a fixed probability of containing a detectable blob, 
which is labeled the detectability fraction /dE3- We 
choose /d such that the halo mass function from the 
simulation reproduces the observed blob number density 
from £14.3.11 Once the halo mass and detectability frac- 
tion are fixed, we can predict the clustering of such halos, 
i.e., the field-to-field variation, directly from the simu- 
lation using counts-in-cells (c-in-c) methodology. Then 
we compare this prediction to the observed variation in 
blob counts over the four survey fields to quantify the 
likelihood that the selected halo mass and detectability 
fraction reproduce the observed blob statistics in Table 

12 Note that /d is a statement about the detectability of the en- 
semble of halos and is not the classically-defined duty-cycle, which 
is simply the fraction of time a blob is on. 
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To link the number density and its variance to the DM 
halo mass, we employ a simple counts-in-cells analysis 
using the halo catalog at z — 2 derived from the ABA- 
CUS N-body code (Metchnik & Pinto, in prep.). This 
simulation has a cubic volume of XlrT 1 co- moving Gpc 
on a side and 1024 3 dark matter particles with raDM = 
l.lxlO 11 M . We adopt the cosmological parameters: 
H = 0.701, n s = 0.96, ^ M = 0.279, Q b = 0.0456 and 
(78 = 0.817. Dark halos are defined using a friends-of- 
friends algorithm with linking length b = 0.16 in units of 
the initial particle spacing. The smallest halos used in 
our analysis consist of 48 particles. Due its large size, this 
simulation is finely tuned to our problem, which requires 
sampling many "cells," ~50 comoving Mpc boxes that 
are roughly the same size and geometry as our survey of 
each of the four fields. 

To constrain the DM halo mass of the blobs, we first 
consider the observed bright/large blob number den- 
sity N = 1.2 per 0.25 deg 2 field, or n = 1.0xl0~ 5 
Mpc" 3 (1.5 blobs in a 50 Mpc box). If all halos con- 
tain a detectable blob (i.e., the detectability fraction /d 
is 100%), then this number density requires halos with 
more than 150 DM particles or 1.7xlO 13 M . We de- 
rive the counts-in-cells distribution of the simulated blobs 
by counting the number of blobs within 10000 randomly 
placed 50 Mpc boxes and b y assuming; a simple halo occu- 
pancy distribution (e.g.. iBerlind fc Weinbe rg 2002) with 
N g (M ^ M min ) = 1 + (M/Mi) a , where M min = 150 DM 
particles or 1.7xlO 13 M , M 1 = 1000 DM particles or 
1.1 x 10 14 Mq, and a — 1 as a fiducial value. This cell size 
is similar to our survey volume, 48.7x48.7x46.8 Mpc, 
for the 0.25 deg 2 FOV with the NB403 narrowband fil- 
ter. While the four survey fields have slightly different 
survey dimensions, the choice of a 50 Mpc box does not 
affect our conclusions. 

Figure [11] shows the simulated bright/large blob 
counts-in-cells distribution with an average of 1.46 blobs 
per cell. The overlayed line represents a Poisson distribu- 
tion with the same average, a reference case in which the 
DM halos containing blobs are not clustered. Thus, com- 
paring the bright/large and Poisson distributions tests 
the null hypothesis that the halos and blobs are not cor- 
related. 

As expected from the large field-to-field variation 
( §4.3.ip . the low and high tails of the simulated distri- 
bution exceed the Poisson counts. Using the Eq. (j2j), we 
derive the variance of this distribution, a v = 1.04, which 
is consistent with our estimate (~1.4a lower limit) from 
the previous section ( §4.3.1)) . This variance corresponds 
to a bias of b — <7bi b / ^dm °° 7 given the rms fluctuation 
of mass ((7dm = 0.15) within a sphere of 31 Mpc radius 
with the same volume as our survey box. 

Is this simulated distribution consistent with the ob- 
served field-to-field variation in n? For /d = 100% 
(shaded histogram), the probability of finding no blobs in 
a cell is as high as41.5%(~2x greater than the Poisson 
probability), while the probability of finding 6 or more 
blobs is 4.4% (~ 11 x greater than Poisson). The poste- 
rior probability of finding at least six blobs in one survey 
field with non-detections in three other fields (the case for 
the bright/large blobs in our sample) is 1.3% (4 x0.415 3 
x 0.044). This posterior probability is ~65 times larger 
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Number of Blobs in the Cell 

Fig. 11. — Counts-in-cells distribution of blobs drawn from the 
ABACUS 1 /i -1 Gpc N-body simulation (Metchnik & Pinto, in 
prep.) (shaded histogram) that have the same average number 
density as our bright, large blob sample (n = l.OxlO -5 Mpc -3 ). 
We count the number of blobs within 10000 randomly placed 
50 Mpc boxes assuming that all massive halos of at least M m j n = 
1.7xl0 13 Mq have detectable blobs (i.e., the detectability fraction 
is 100%). The simulated distribution has a Poisson-noi se corr ected 
variance o~ v = 1.04, consistent with our estimates from £14.3.11 The 
thick solid histogram represents the distribution for a detectability 
fraction of 12.5% and a blob halo mass of M m j n = 0.5xl0 13 Mq. 
The observed number density and the strong field-to-field varia- 
tion of the bright, large blobs at z = 2.3 suggest that they reside 
in halos that are likely to evolve into the ~ 10 14 Mq halos typical 
of present-day clusters of galaxies. The dots represent the Poisson 
noise with a mean of 1.46. The inset shows the same distribution 
but for a logarithmic scale. 

than 0.02% for the Poisson distribution. Therefore, we 
reject the null hypothesis that the dark matter halos and 
blobs are uncorrelated. 

Although the posterior probability is much higher than 
Poisson, is it statistically acceptable as a halo mass con- 
straint? To test if we can reject the posterior probability 
of 1.3%, we consider an extreme case that maximizes our 
posterior statistics: P = 4p 3 (l— p), the probability of not 
finding any blob in three survey fields, but finding any 
number of blobs in the one remaining field. Here, p rep- 
resents the probability of finding zero blobs in one survey 
field, and P has a maximum value of 42.2% when p = 
0.75. Therefore, we should compare our posterior prob- 
ability of 1.3% with this extreme case, not with 100%. 
The posterior probability is only ~34x smaller than this 
maximum probability. Therefore, we conclude that our 
measurements are consistent with the halo model assum- 
ing M halo ~10 13 M and f D = 100%. 

We have now established that the observed strong field- 
to-field variation in blob counts is not surprising if mas- 
sive dark matter halos with Mh a io ^ 10 13 M© always pro- 
duce detectable Lya blobs. Here we investigate whether 
it is possible to obtain consistency with lower values of 
halo mass by changing the blob detectability fraction. 
We derived the maximal halo mass for the observed blob 
number density by assuming that all halos more massive 
than Mhaio °° 1.7xlO 13 M have detectable blobs. If we 



lower /d to 50%, 25%, and 12.5%, while increasing the 
number density of halos capable of hosting blobs by 2, 4, 
and 8 times to keep the abundance of the observed blobs 
constant, the threshold halo mass M m i n decreases to 1.2, 
0.8, and 0.5 x 10 13 M , respectively. At lower halo mass, 
we naturally obtain a weaker field-to-field variation and a 
less prominent high-end tail in the counts-in-cells distri- 
bution (i.e., at N > 6) than for /b = 100%. To illustrate 
this trend, we also show the counts-in-cells distributions 
for / D = 12.5% (thick solid line) in Figure HU Table 
[71 summarizes the counts-in-cells statistics, the field-to- 
field variance (cr v ), and the posterior probability of ob- 
serving 6 or more blobs in only one of the survey fields for 
the different /b values. For comparison, we also list the 
counts-in-cells statistics and the posterior probabilities 
for the Poisson distributions. 

Higher detectability fractions (and more massive ha- 
los) produce a field-to-field variance about the observed n 
that is more consistent with the observed variance. Halo 
models with > 10 13 M halos and /b > 50% work best. 
Lower detectability fractions (e.g., ~ 12%) require some- 
what lower halo masses (~ 5 x 10 12 Mq) to reproduce n, 
but the resulting variance in n is lower and further from 
the observed value. However, the effects of lowering the 
detectability fraction and halo mass are not large enough 
to put strict lower limits on the halo mass: for the low- 
est /d considered, 12.5%, we predict a v = 0.76 (76%) 
and a posterior probability of 0.45%, which is only ~ 3x 
lower than in the /b = 100% case. It is possible, but less 
likely, that the blobs occupy a few x 10 12 M halos if /b 
is much lower than 12.5%. 

A detectability fraction of /b = 12.5% implies a short 
blob lifetime, only r < 350 Myrs at z = 2 — 3. In princi- 
ple, we could adopt still lower values (/b <C 10%), down 
to the limit where the Lya blobs live only a few tens 
of Myrs, and thus lower halo masses (<1O 12 M ), in or- 
der to find the point at which we can reject the assumed 
halo mass. However, the limited mass resolution of our 
N-body code does not allow us to resolve smaller halos. 
Therefore, to put tighter constraints on the halo mass 
requires improving the blob statistics by extending sur- 
veys to larger volumes and/or creating higher resolution 
simulations. 

For now, we conclude that bright/large Lya blobs 
are most likely to reside in massive dark halos with > 
10 13 M that have detectable blobs more than ~ 50% 
of the time. Interestingly, these halo mass estimates 
agree with the dynamical masses, M^ yn = 10 12 -10 13 M , 
derived from the wid th of Lya lines in similar blobs 
(jMatsuda et al.l 12006). However, special caution is re- 
quired in using the Lya line width as a mass proxy, be- 
cause the radiative transfer effects on the line width are 
poorly understood. 

For the entire blob sample (with counts of 14, 5, 3, 1 
for the four survey fields), the required halo mass for /b 
= 100% is M min = 0.8 x 10 13 M Q . The fractional variance 
is <j v = 0.76, again consistent with the observed value. 
Because the halo mass function is steep at the high mass 
end, M m i n for the entire sample is similar to that derived 
for the bright/large blobs alone. Therefore, the resulting 
halo mass is insensitive to the blob selection cuts, and 
the statistics for our entire sample are consistent with 
blobs occupying halos of ~ 1O 13 M . 
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Halos of ~ 10 13 Mq lie in the high mass tail of the 
halo mass distribution at z = 2.3. How massive will 
these halos be today? Because our simulations were not 
ru n beyond z = 1, w e consider the N-body simulation 
of iMaccio et al.l ([2008[ h which extends to z = 0. Their 
Miaio ^ 10 13 M Q halos grow in mass by 2-10 x (with 
an average factor of 5.2 ± 2.4 increase) from z = 2.3 to 
now. Therefore, it is likely that blobs are sites for the 
formation of the brightest galaxies in what will become 
the typical halos of rich clusters (Mh a io °° 1O 14 M ) at 
the present epoch. 

4.4. Large-Scale Environment of Blobs 

The clustering of Lya blobs in the CDFS and the in- 
ferred large mass of their individual halos imply that 
blobs inhabit over dense regions. Here we test this hy- 
pothesis using the much larger population of compact 
Lya emitters (LAEs) in the CDFS to trace large-scale 
structure over tens of comoving Mpc and thus character- 
i ze th e blob environment independently. Mats uda et al.l 
(2005) show that their blobs lie near the intersection of 
large -scale filamentary stru ctures in the SSA22 overden- 
sity. Pre scott et al.l ([2008) use the surface number den- 
sity of LAEs to show that a gia nt blob, originall y identi- 
fied via its strong MIR emission (|Dev et al.ll2005l ), resides 
in a region ~ 3 x more dense than the edge of their survey 
field. We apply counts-in-cell methodology to the LAE 
spatial distribution in the CDFS to quantify the scale 
over which any structure is coherent. We then identify 
over-densities of LAEs relative to their average number 
density in the field and compare them to the spatial dis- 
tribution of blobiQ 

We select a sample of LAEs in the CDFS (TV ~ 200 
with NB < 25.0, excluding the blobs) following the first 
step in our blob selection procedure (^3j), but using a 2"- 
diameter aperture to maximize the S/N of fainter point 
sources. In Table [TJ we list the 5a limiting magnitudes, 
which are determined by measuring fluxes within the ran- 
domly placed apertures in the sky background region. 

Using this LAE catalog, we test whether there is struc- 
ture over various spatial scales. For each scale, we count 
the number of LAEs within circular cells of that scale 
radius randomly placed over the field. Then we compare 
this counts-in-cells distribution with a Poisson distribu- 
tion using Kolmogorov-Smirnov (K-S) statistics to check 
if this distribution deviates from Poisson noise. We re- 
peat this test with different cell sizes ranging from 1' to 
8'. The counts-in-cells (c-in-c) distributions deviate sig- 
nificantly (at > 95.0% confidence) from Poisson noise for 
cells larger than 5', indicating coherent large-scale struc- 
ture over scales of at least ~ 8 comoving Mpc. 

Figure [[2] shows the LAE distribution in the CDFS 
field. To estimate the over density of LAEs, we overlay 
the surface over density = — £)/£ contours on the 
map. Here E represents the average surface density esti- 
mated from all the LAE candidates over the entire field. 
The overdensity maps are smoothed using the adaptive 
kernel smoothing method developed by Ebe ling et al.l 
(|2006[ ). The FWHMs of the Gaussian filters adopted by 

14 We do not detect any significant large-scale structure in the 
other three survey fields. Because survey depth and seeing vary 
over the fields, it is not clear whether the non-detections arise from 
a real absence of structure or poorer sensitivity. 
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Fig. 12. Surface density map (~30'x30'= 50Mpcx50Mpc) of 
compact Lya emitters in the CDFS. The small and big dots repre- 
sent compact and extended Lya: emitters (i.e., LAEs and LABs), 
respectively. The contours are spaced for 5-£ = (E — = —0.7, 

—0.5, 0.0, 0.5, 1.0, where S is the the average LAE surface density 
of the whole region. The Lya blobs are located in the overdense re- 
gion traced b y the LAEs. The three crosses represent overdensities 
identified by Salimbeni et al. (2009) using photometric redshifts 
(zphot = 2.23,2.28,2.28). A counts-in-cells analysis of the LAEs 
confirms coherent large-scale structure over scales of at least ~ 8 
comoving Mpc. 

this algorithm ranges from 2.8' to 5.2' depending on the 
local surface density of LAEs within in the field. We iden- 
tify a belt-like large-scale structure elongated from north 
to south. The blob candidates are preferentially located 
within or at the boundary of this structure. Therefore, 
we conclude that the Lya blobs in the CDFS indeed trace 
overdensities i n the early unive rse. 

Recently, I Salimbeni et al.l (|2QQ9h systematically 
searched for galaxy overdensities within the GOODS- 
South field up to z ~ 2 using photometric redshifts that 
have a rms accuracy of Az/(1 + z) = 0.03. In Figure [T2l 
we show their three overdensities (crosses) with redshifts 
similar to our survey (z p hot — 2.23-2.28). While the 
GOODS-South field does not include the densest part of 
our surface density map, their overdensities each include 
19-23 members and are located near the boundary of 
the LAE overdensity. Spectroscopic confirmation of 
this system is required to relate these systems with the 
structures revealed by our Lya blobs and emitters. 

Further studies are required to establish that the CDFS 
overdensity a nd similar overdensiti es observed by other 
authors (e.g., iPrescott et al.l 120081 ) are in fact "proto- 
clusters." For now, we note that structure in the CDFS 
at z = 2.3 is coherent over at least 8 comoving Mpc and 
is likely to contain tens of ~ 10 13 M blob halos, or at 
least ~ 10 14 Mq worth of mass. Were this overdensity to 
grow typically and virialize by z = 0, then its halo mass 
today would be like that of the rich Coma cluster, i.e., ~ 
10 15 M«. 
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5. CONCLUSIONS AND SUMMARY 

To understand what Lya blobs will become in the 
present-day universe requires that we first constrain their 
number density, clustering, and large-scale environment. 
In order to obtain unbiased measures of these quanti- 
ties, we target four ^jS0 f x30 f field s, Chandra Deep Field 
South (CDFS; Brandt et al. 2001) Chandra Deep Field 
North (CDFN ; Giacconi et al. 200 2]), and two COSMOS 
(|Scoville et all 120071 : Koekem oer et all I2007D subfields, 
with the NOAO Mayall and Blanco-4m telescopes and 
a custom narrowband filter designed for Lya at z ~ 2.3. 
The total area covered by our survey is 1.2 deg 2 . Our 
sensitivity and selection criteria a re comparable to tha t 
of the largest previous blob survey (Ma tsuda et al.l l2Q04). 

We discover 25 Lya blobs with Lya luminosities of 
Lh y a= 0.7-8 x 10 43 ergss -1 and isophotal areas of A{ so 
= 10-60 n". The transition from compact Lya emitters 
(LAEs; 

°° a few n'O to the extended Lya blobs (A[ so 
> 10D") is continuous, suggesting that these two types 
of sources are not distinct and that whatever mechanism 
or mechanisms power Lya blobs work over a wide range 
of luminosity and spatial extent. 

Surprisingly, we find the majority of blobs (16/25) in 
one survey field, the CDFS. The six brightest (L\, y0i > 
1.5 xlO 43 ergss -1 ) and largest (A- lso > 16D") blobs are 
discovered only in CDFS, indicating a strong field-to- 
field variation. Using a simple analytic approximation for 
the underlying fluctuations of the blob number density, 
we find that these large/bright blobs have a field-to-field 
variance of a v > 1.5 (150%) about their number density 



of 



1.0 



-0.6 ' 



10 Mpc . This variance is large, 



significantly higher than that of unresolved Lya emitters 
(a v -0.3 or 30%). 

To constrain the mass of the dark matter halo around 
each Lya blob, we compare the number density and clus- 
tering of blobs with the counts-in-cells distribution of 
halos predicted from a 1 h~ 1 Gpc cosmological N-body 
simulation. At z = 2.3, n implies that bright, large blobs 
could occupy halos of Mh a io ^ 10 13 M if most halos 
contain a detectable blob, i.e., the detectability fraction 



is > 50%. The predicted variance in n is consistent with 
that observed and corresponds to a bias of 7. Lower 
detectability fractions (e.g., ~ 10%) require somewhat 
lower halo masses (~ 5 x 1O 12 M ) to reproduce n, but 
the resulting variance is lower and further from the ob- 
served value. Blob halos lie at the high end of the halo 
mass distribution at z = 2.3 and are likely to evolve into 
the ~ 10 14 Mq halos typical of galaxy clusters today. 

The clustering and inferred halo mass of blobs suggest 
that they lie in over dense environments. The spatial dis- 
tribution of LAEs confirms this hypothesis: a counts- 
in-cells analysis of the CDFS reveals coherent large-scale 
structure over scales of at least ~ 8 comoving Mpc where 
both the LAEs and blobs cluster. 

Given the strong field- to-field variance of Lya blobs, 
one must be cautious in comparing blob number den- 
sities and luminosity functions among different surveys. 
We construct a reliable luminosity function of Lya blobs 
from a deep, blind narrowband survey. Larger volume 
blob surveys, combined with large volume and/or higher 
resolution N-body simulations, will improve the con- 
straints on blob halo mass and detectability fraction, 
thus discriminating among the possible energy sources 
of the extended Lya emission. 
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TABLE 1 

Narrowband Observations 



Field 


R.A. 


Dec. 


Observing Date 


Telescope 


Exposure 


Depth a 


Survey Area 


Seeing 


Pixel Scale b 




(J2000) 


(J2000) 






(hour) 


(AB mag) 








CDF-S 


03:32:27.8 


-27:47:56 


Nov. 2007 


Blanco 4m 


10 


25.65 


31.'6 x 31.'6 


1.0" 


0"27 


CDF-N 


12:36:50.9 


62:11:48 


May 2007 


Mayall 4m 


10 


25.27 


29.'5 x 29.'5 


1.3" 


0'/30 


COSMOS1 


09:59:16.9 


01:55:19 


Feb. 2009 


Blanco 4m 


7.7 


25.27 


36.'0 x 36.'3 


1.2" 


0"27 


COSMOS2 


09:59:16.8 


02:31:19 


Feb. 2009 


Blanco 4m 


7.2 


25.25 


36.'0 x 36.'3 


1.2" 


0"27 



a 5<t detection limit for 2 // -diameter aperture. 

k Pixel scales of final combined images, which are determined by the largest pixel scale between narrowband and broadband images. 



TABLE 2 
Summary of Broadband Images 



Field 


Band 


Effective Wavelength 


Band Width 


Depth 


Seeing 


Telescope and Instrument 


Reference 






(A) 


(A) 


(AB mag) 


(arcsec) 






CDF-S 


U 


3505 


626 


26.0 a 


1.05 


ESO 2.2m WFI 


Gawiser et al. (2006b) 




B 


4600 


915 


26.9 


0.95 


ESO 2.2m WFI 




CDF-N 


U 


3648 


387 


27.1 b 


1.26 


KPNO MOSAIC I 


CaDak et al. (2004) 




B 


4428 


622 


26.9 


0.71 


Subaru Suprime-Cam 


COSMOS 


U 


3798 


720 


26. 4 b 


0.90 


CFHT Megaprime 


CaDak et al. (2007) 




B 


4460 


897 


27.3 


0.95 


Subaru Suprime-Cam 



a CDF-S: 5cr detection limit corrected for infinite aperture. 

b CDF-N and COSMOS 5a detection limit for Skiameter aperture. 



TABLE 3 

Properties of Lya Blob Candidates 



ID 


R.A. 
(J2000) 


Dec. Ltot L iso 
(J2000) (lO^ergss- 1 ) (10 43 ergs s" 1 ) 


Area 
(□") 


Note 


Extended Chandra Deep Field South 



CDFS-LAB01 


03 


32 


36.1 


-28 00 


54.5 


7.81 


± 


0.28 


8.00 


61.6 


± 


3.7 


confirmed 


CDFS-LAB02 


03 


33 


20.6 


-27 41 


10.8 


3.26 


± 


0.17 


2.88 


37.9 


± 


3.3 


confirmed 


CDFS-LAB03 


03 


31 


52.1 


-27 54 


54.6 


3.23 


± 


0.29 


2.72 


43.2 


± 


4.4 




CDFS-LAB04 


03 


32 


05.9 


-27 37 


16.7 


2.64 


± 


0.15 


2.40 


30.9 


± 


2.1 


confirmed 


CDFS-LAB05 


03 


31 


48.7 


-27 52 


23.3 


1.72 


± 


0.39 


1.14 


21.3 


± 


6.2 




CDFS-LAB06 


03 


32 


19.6 


-27 47 


30.8 


1.53 


± 


0.11 


1.57 


16.5 


± 


1.8 




CDFS-LAB07 


03 


32 


03.2 


-27 45 


25.0 


1.48 


± 


0.17 


1.56 


18.5 


± 


2.5 




CDFS-LAB08 


03 


31 


43.7 


-27 38 


32.9 


1.21 


± 


0.22 


1.32 


19.3 


± 


4.5 




CDFS-LAB09 


03 


32 


44.5 


-27 43 


10.2 


1.17 


± 


0.15 


0.87 


14.3 


± 


2.0 




CDFS-LAB10 


03 


32 


37.4 


-28 02 


05.7 


1.14 


± 


0.25 


0.71 


12.6 


± 


2.6 


confirmed 


CDFS-LAB11 


03 


32 


43.2 


-27 42 


58.3 


1.09 


± 


0.08 


1.02 


10.3 


± 


1.3 




CDFS-LAB12 


03 


32 


06.7 


-27 44 


55.2 


1.08 


± 


0.16 


0.73 


13.3 


± 


2.4 




CDFS-LAB13 


03 


32 


32.7 


-27 39 


06.3 


1.03 


± 


0.09 


0.94 


12.0 


± 


1.5 




CDFS-LAB14 


03 


32 


32.2 


-27 41 


27.2 


0.98 


± 


0.09 


0.93 


12.8 


± 


1.5 


confirmed 


CDFS-LAB15 


03 


31 


38.6 


-27 43 


07.0 


0.96 


± 


0.08 


0.81 


11.8 


± 


1.1 




CDFS-LAB16 


03 


33 


05.8 


-27 57 


40.0 


0.69 


± 


0.18 


0.62 


10.3 


± 


3.0 





Chanda Deep Field North 



CDFN-LAB01 


12 35 35.2 


+62 14 28.5 


1.31 ± 0.11 


1.30 


13.2 ± 1.8 






CDFN-LAB02 


12 35 30.2 


+62 01 39.2 


1.15 ± 0.12 


0.94 


11.7 ± 1.7 


mar^ 


pnal? 


CDFN-LAB03 


12 36 09.9 


+61 57 16.6 


1.01 ± 0.12 


0.94 


12.2 ± 1.9 


mar^ 


pnal? 


CDFN-LAB04 


12 36 59.2 


+62 24 35.2 


0.99 ± 0.15 


0.80 


13.0 ± 2.0 






CDFN-LAB05 


12 38 18.4 


+62 04 04.1 


0.72 ± 0.13 


0.61 


11.1 ± 1.9 










Cosmic Origins 


Evolution Survey Field 








COSMOS-LAB01 


09 59 23.9 


+01 55 11.7 


1.24 ± 0.19 


0.87 


14.6 ± 5.2 


marj 


^inal? 


COSMOS-LAB02 


09 59 14.2 


+01 48 43.1 


0.92 ± 0.19 


1.05 


11.9 ± 2.3 


mar^ 


pnal? 


COSMOS-LAB03 


09 58 12.8 


+01 52 56.2 


0.84 ± 0.16 


0.70 


12.5 ± 2.5 






COSMOS-LAB04 


09 58 49.6 


+02 30 49.9 


1.01 ± 0.13 


0.92 


11.8 ± 2.1 







Note. — The note column indicates whether the redshifts of the candidates are spectroscopically confirmed 
with the Lyo; and/or Ha line and whether they are marginally resolved (labelled "marginal"). 
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TABLE 4 
Blob Statistics 





CDFS 


CDFN 


COSMOS1 


COSMOS2 


A iso > 16 □" 


8 











A iso = 10-16 □" 


8 


5 


3 


1 


^blob (bright /large sample) a 


6 











^blob (entire sample) 


14 


5 


3 


1 


Effective area [deg 2 ] 


0.270 


0.241 


0.352 


0.362 



a Nhiob is the number of blobs corrected to the same depth and seeing for all 
four survey fields. We use iVbiob for estimating the number density, its variance, 
and the blob halo mass in ^4.3.11 



TABLE 5 

FlELD-TO-FlELD VARIANCE AND NUMBER DENSITY OF BLOBS 





Variance o v Number per 0.25 deg 2 N 




Min Best Max Min Best 


Max 


Bright, Large Blobs with A iso > 16 


Joint 
Marginal 


1.49 (0.58) 0.27 

3.11 (0.99) 0.40 (0.17) 1.15 


3.27 (8.30) 


Entire Sample with A iso > 10 □" 


Joint 
Marginal 


0.26 (0.10) 0.87 5.21 2.27 (1.56) 4.76 
0.57 (0.31) 1.20 3.51 2.91 (2.13) 4.76 


21.7 

10.3 (27.1) 



Note. — The number density N represents the number of blobs per 0.25 
deg 2 . In the parentheses, we list the 2a limits. 



TABLE 6 

CUMULATIVE LUMINOSITY FUNCTION OF LYq; BLOBS 



Log(L Lya ) n(> L) [xlO- 5 Mpc- 3 ] 





CDFS a 


All b 


SSA22 C 


43.90 






1.54 ± 1.09 


43.70 


0.83 ±0.83 


0.18 ±0.18 


2.31 ± 1.33 


43.50 


2.50 ± 1.44 


0.55 ±0.32 


3.08 ± 1.54 


43.30 


3.33 ±1.67 


0.73 ±0.37 


5.38 ±2.04 


43.10 


5.83 ±2.20 


1.47 ±0.52 


10.77 ±2.88 


42.90 


12.49 ±3.22 


4.23 ±0.88 


18.46 ±3.77 


42.70 


13.32 ±3.33 


4.59 ±0.92 


26.92 ±4.55 



a LF from blobs in CDFS. 

b LF from all four survey fields (CDFS, CDFN, COSMOS1, and 

COSMOS2). 

c We also list the LF from[Matsuda et al. (2004) for comparison. 
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TABLE 7 

COUNTS-IN-CELLS DISTRIBUTION FOR BRIGHT, LARGE BLOBS 



N-Body Simulation Poisson Distribution 



fD 


M min (M ) 


a v 


P(N = 0) 


P(N ^ 6) 


Probability 


P(N = 0) 


P(N ^ 6) 


Probability 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


100% 


1.66X10 13 


1.03 


41.5% 


4.40% 


1.25% 


23.3% 


0.440% 


0.022% 


50% 


1.16X10 13 


0.93 


38.5% 


3.65% 


0.83% 


22.9% 


0.370% 


0.018% 


25% 


7.97X10 12 


0.83 


35.8% 


3.27% 


0.60% 


23.4% 


0.300% 


0.015% 


12.5% 


5.32X10 12 


0.76 


35.6% 


2.51% 


0.45% 


23.8% 


0.370% 


0.020% 



Note. — (1) detectability fraction, (2) dark matter halo mass, (3) variance [Eq. J2J], (4) probability of non- 
detection, (5) probability of finding 6 or more blobs, (6) posterior probability of finding 6 or more blobs in one field 
and none in the three other fields, (7) (9) same as (4) (6) but for Poisson distribution. 



100 



Q 

O 




Isophotal Area of 35 SSA22 Blobs [arcsec 2 ] 




10 4a 10* 
Luminosity of 35 SSA22 Blobs [ergs s" 2 ] 



Fig. 13. — Comparison in L\, ya and A- lso between our 16 CDFS blobs (horizontal bars) and the 35 Matsuda et al. (2004) blobs in the 
SSA 22 field (filled circle s). The vertical and horizontal dashed lines indicate the size selection criteria for the extended sources adopted 
by [Matsuda et al. ( 2004) and by our work, respectively. The dotted lines re present the average ra tios (~0.61 and 1.0) between originally 
reported and the recovered values in A- lso and L^ ya . The recovery test of the Matsuda et al. (2004) blobs when pasted into our CDFS field 
shows that we can detect all blobs like theirs in our survey. 



APPENDIX 

COMPARISON WITH MATSUDA ET AL. Lva BLOB SAMPLE 

It is difficult to compare the properties of Lya blob samples among different survey s because of non - unifo rm selection 
criteria and different imaging depth. Here we compare our z = 2.3 sampl e to tha t of Matsud a et al.l (j2004f ) at z = 3.1. 
We repeat t he procedures for the recovery test described in $3] (see also lYang et al.ll2009[ ), pasting thumbnail images 
provided bv lMatsuda et al.l ((2004) into our CDFS image. Because the narrowband filter bandwidths of the two surveys 
are similar, we do not make any correction for the difference in filter transmission, but we scale the apparent size and 
the surface brightness accounting for the different redshifts. 

Figure [13] sho ws the Lya luminosi ty and isophotal area recovered from this test as a function o f the input Z/T jV n , and 
A[ so for the 35 lMatsuda et all (j2QQ4f ) blobs. Due to the slightly shallower depth of our survey, the Ma tsuda et al.l blobs 
would look smaller by a factor of 61% than the originally reported sizes. Note that one expects a 86% decrease in area 
purely from the differences in angular diameter distance between two survey reds hifts. The vertical an d horizontal 
dashed lines indicate the size selection cr i teria f or the extended sources adopted by Matsud a et al.l (|2QQ4f ) and by our 
work, respectively. Most Matsu da et al.l (j2004f ) blobs are recovered as larger than 10 confirming the capability 
of our survey for detecting them. The line luminosities are also recovered well, so there is no bias in the luminosity 
measurements. 



